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en ■ Abstract: We provide the first theoretical study of a novel variable, aT, proposed in 

. ref. [1] as a more accurate probe of the region of low transverse momentum for the Z 

. boson pt distribution at hadron colliders. The ax is the component of px transverse to 

0^ i a suitably defined axis. Our study involves resummation of large logarithms in a-r up to 

' the next-to-leading logarithmic accuracy and we compare the results to those for the well- 
known Pt distribution, identifying important physical differences between the two cases. 

^ ' We also test our resummed result at the two-loop level by comparing its expansion to order 

■ a with the corresponding fixed-order results and find agreement with our expectations. 
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1. Introduction 

The study of W and Z boson production at hadron colhders via the Dreh-Yan process [2] 
has formed a very significant part of particle phenomenology over a period of several 
years [3~5] . In the era of the LHC these studies continue to occupy an important role for 
a variety of reasons. For instance an accurate understanding of the production rates and 
transverse momentum {px) distributions of the W and Z, or equivalently those of lepton 
pairs obtained from gauge boson decays, can be used for diverse purposes. These range 
from more prosaic applications such as luminosity monitoring at the LHC to measurement 
of the W mass and perhaps most interestingly for discovery of new physics, which may 
manifest itself via the decay of new gauge bosons to lepton pairs. 

In this paper we shall concentrate on the pT spectrum of the Z boson and related 
quantities. The Z boson pT distribution has in fact received considerable theoretical and 
experimental attention in the past but there remain aspects where it is desirable to have 
an improved understanding of certain physical issues. One such important issue is the role 
of the non-perturbative or "intrinsic" component which has a sizable effect on the pT 
spectrum at low px (see e.g. refs. [6-8]). The fact that the incoming quarks/anti-quarks 
which fuse to form gauge bosons are part of extended objects (protons or anti-protons) and 
have interactions with other constituents thereof generates a small transverse momentum 
(that can be viewed as the Fermi motion of partons inside the proton) and which a priori 
one might expect to be of order of the QCD scale Aqcd- 
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Since the intrinsic kx has a non-perturbative origin it cannot be computed within 
conventional methods of perturbative QCD. One can however model the intrinsic fey as an 
essentially Gaussian smearing of perturbatively calculated px spectra and hope to constrain 
the parameters of the Gaussian by fitting the theoretical prediction to experimental data. 
An example of this procedure is provided by the work of Brock, Landry, Nadolsky and 
Yuan (BLNY) who proposed a non-perturbative Gaussian form factor that in conjunction 
with perturbative calculations was able to describe both Tevatron Run-1 Z data as well 
as lower Drell-Yan data [8]. Alternatively one may use a Monte Carlo event generator 
such as HERWIG++ [9] to phenomenologically investigate the same issue. As discussed 
in ref. [10] these studies yield somewhat larger than expected values for the mean intrinsic 
kx per parton and also reveal a dependence of this quantity on the collider energy which 
features are desirable to understand better. 

Additionally, as pointed out by Bcrgc et al. [11], a phenomenological study in semi- 
inclusive DIS processes (SIDIS) for small Bjorken-x, x < 10~^, suggests an effective depen- 
dence of the non-perturbative form factor on x.^ Extrapolating the effect to the LHC where 
such small-x values become relevant, one may expect to see significantly broader Higgs and 
vector boson px spectra than one would in the absence of small-x effects. Tevatron studies 
with high rapidity vector boson samples may help to provide further information on the 
role, if any, of small- .t broadening. 

Given the importance of the studies we have mentioned above, in the context of the 
LHC and the precise determination of px spectra there, it is important to have as thorough 
a probe of the low px region of Z boson production as is possible. Investigations carried out 
using the conventional Z boson px spectrum mainly suffer from large uncertainties arising 
from experimental systematics, dominated by resolution unfolding and the dependence on 
Px of event selection efficiencies, as discussed in detail in ref. [1]. 

An interesting observable that is less sensitive to the above effects and also essentially 
insensitive to the momentum resolution of leptons produced by the decaying vector boson 
was proposed in ref. [1]. This variable is just the transverse component, ax, of the px with 
respect to the lepton thrust axis, which we shall define more precisely later. It has been 
suggested that the ax variable will be experimentally better determined at low ax than 
the standard px variable and hence it would make a more accurate probe for issues such as 
our understanding of initial state radiation and the precise role of intrinsic kx, including 
issues such as potential small-a; broadening. 

Before one can access information on non-perturbative effects however, it is of vital 
importance to have a sound perturbative estimate of the observable at hand. In the low 
Px region of interest to us one is dealing with the emission of soft and/or collinear glu- 
ons which is logarithmically enhanced. The resummation of large logarithms of the form 
1/px (a" In*" M'^ /p\) , where M is the lepton-pair invariant mass and m < 2n— 1, has been 
a subject of interest over decades [13-18] and has now been carried out to next-to-next-to 
leading logarithmic (NNLL) accuracy [19]. After matching such resummations with fixed- 

^This X dependence may be merely an eflfective parametrisation of missing perturbative BFKL effects. 
Another observable, the vector pr of the current hemisphere in the DIS Breit frame which can be used to 
investigate this x dependence, was suggested in ref. [12]. 
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order estimates to NLO accuracy one has a state-of-the art theoretical prediction for the 
perturbative region. 

In the present paper we address the issue of resummation for the ax variable. We point 
out that there are some similarities to resummation for the pr distribution but also some 
important differences that manifest themselves in the shape of the resummed distribution. 
We aim to provide a next-to-leading logarithmic (NLL) resummation that we envisage 
could be extended to NNLL level subsequently. The NLL resummed form we provide here 
can however already be used after matching to full next-to-leading order (NLO) results for 
accurate phenomenological studies of ay- 

This paper is organised along the following lines: we begin by discussing the definition 
of the ax and its dependence on multiple soft gluon emission which is important at low 
ut, where one encounters large logarithms. In the following section we sketch a leading- 
order calculation for the ut distribution, consigning details to an appendix, which helps to 
illustrate some features of the qt distribution such as the precise origin of logarithmically 
enhanced terms. In the subsequent section we carry out a resummation of the logarithms 
of qt to NLL accuracy pointing out the relation to our recent work on azimuthal jet 
decorrelations [20]. Next we identify a relationship between the ar and pj^ distributions at 
fixed order and check this relationship with the help of a numerical fixed-order calculation 
using the program MCFM [21], which is a non-trivial test of our resummation. We conclude 
by pointing out the possibilities for further work which involve a possible extension of our 
resummation to NNLL accuracy (as has been done for the px distribution [19]) as well as 
matching to the MCFM results and phenomenological investigation once final experimental 
data becomes available. 



2. Definition of a.T and soft limit kinematics 

We will be concerned in this paper with large logarithms in the perturbative description 
of the ax variable and their resummation. Since these logarithms have their origin in 
multiple soft and/or collinear emissions from the incoming hard partons we need to derive 
the dependence of the ax on such emissions. In this section therefore we define precisely 
the ax and obtain its dependence on the small transverse momenta kt of emissions. 

We recall that we are considering the production of Z bosons via the Drell-Yan (and 
QCD Compton) mechanisms which subsequently decay to a lepton pair. The ax is the 
component of the lepton pair (or equivalently Z boson) px transverse to a suitably defined 
axis. The precise definition of the lepton thrust axis as employed in ref. [1] is provided 
below: 

Pa - Pt2 , ^ 

\Ptl -Pt2\ 

where pn and pt2 are the transverse momenta of the two leptons and thus n is a unit vector 
in the plane transverse to the beam direction. It is straightforward to verify that this is 
the axis with respect to which the two leptons have equal transverse momenta. 

We now consider multiple emissions from the incoming partons which (neglecting the 
intrinsic kx) are back-to-back along the beam direction. From conservation of transverse 
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momentum we thus have pu + pt2 = — j ku which means that the lepton pair or Z boson 
Pt is just minus the vector sum of emitted gluon transverse momenta ku, where we refer 
to the momentum transverse to the beam axis. To obtain the dependence of ot on the ku 
we wish to find the component of this sum normal to the axis defined in eq. (2.1). The 
axis is given by (writing pt2 in terms of pn and ku) 

^ ^ 2pti + ku _ pti_ ^2 2) 

\2pti + Zikti\^ \Pti\' 

where to obtain the last equation we have neglected the dependence of the axis on emissions 
ku- The reason for doing so is that we are projecting the vector sum of the ku along and 
normal to the axis and any term O {ku) in the definition of the axis impacts the projected 
quantity only at the level of terms bilinear or quadratic in the small ku- Such terms can 
be ignored compared to the leading linear terms ~ ku that we shall retain and thus to our 
accuracy the axis is along the lepton direction.^ 

We can parametrise the lepton and gluon momenta in the plane transverse to the beam 
as below: 

m=nii,o), ^^^^^ 

hi = hi (cos 4>i, sin (f)i) , 

where (f)i denotes the angle made by the i^^ emission with respect to the direction of lepton 
1 in the transverse plane. It is thus clear that, expressed in these terms, the transverse 
component of the Z boson px is simply — Y^- hi sin <j)i and one has 



hi sin ( 



(2.4) 



We note immediately that the dependence on soft emissions is identical to the case of 
azimuthal angle Acf) between final state dijets near the back-to-back region A(f) ~ vr, for 
which resummation was carried out in ref. [20]. This is not surprising since the component 
of the Z boson pT, transverse to the axis defined above, is proportional to vr — Acj), where 
A(p is the angle between the leptons in the plane transverse to the beam. The other 
(longitudinal) component of Z boson px, a^, is proportional in the soft limit to pn — Pt2 
the difference in lepton transverse momenta.^ This kinematics is summarised in fig. |l], 
which shows final state momenta in the transverse plane. Together with pn and pt2, the 
two lepton transverse momenta, we have displayed the vector boson transverse momentum 
Pt, the axis h defined in eq. (2.1), and the two transverse momentum components Sl 



and St- From the figure it is also clear that the angle vr — Acp, also indicated, is well 
approximated by | ot I /|pi2 1 ~ clt/ph- In the case of dijet production the kinematics is the 
same, with pn and pt2 representing the transverse momenta of the two highest-pj jets. 

^To be more precise the recoil of the axis against soft emissions, if retained, corrects our result only by 
terms that vanish as — >■ 0. Such terms are beyond the scope of NLL resummation but will be included 
up to NLO due to the matching. 

^For the case of dijet production this pt imbalance has also been addressed via resummation in ref. [22] 
which to our knowledge is the first extension of the pr resummation formalism to observables involving 
final state jets. 
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Figure 1: A di-lepton event in tlie plane transverse to tlie beam. Besides final state momenta, 
the two components of vector boson transverse momentum olt and St are displayed, as well as 
TT — A0, where A0 is the angle between the transverse momenta of the two leptons. See text for a 
full description. 



Since it is possible to measure more accurately the lepton angular separation com- 
pared to their pt imbalance (where momentum resolution is an issue) , one can obtain more 
accurate measurements of ot as compared to ai or the Z boson px which is given by 
aj< + a\ [1]. The resummation that we carry out here will be similar in several details to 
those of refs. [20, 22] but simpler since the final state hard particles are colourless leptons. 
We carry out a leading order calculation of the observable in question in the following 
section and finally the resummation. 



3. Leading order result 

Here we shall mention how to compute the logarithmically enhanced terms in ax at leading 
order in a^, with details of the derivation left to the appendix We shall highlight the 
origin of the double and single logarithms and relate them to the corresponding logarithms 
in the standard pT distribution at the same order. The discussion here should facilitate 
an understanding of the resummation we carry out in the next section and the results of 
subsequent sections. 

At Born level we have to consider the process Pi+ P2 = h + h where pi, p2 and li, I2 
are the four momenta of incoming partons and outgoing leptons respectively. The squared 
matrix element for the Drell-Yan process for lepton pair production via Z decay is [23] 

A^dy(^i, ^2) = ]| ^ («' ^w, M\Ml) [Ai A, {tl + tl) + Bi B, (t? - tj)] , (3.1) 
with the precise forms of the electroweak constant coefficients Ai,Ag and Bi,Bq as well 



as Q reported in eq. ( B.3 ). Henceforth we shall suppress the dependence of Q, which has 
dimension M~^, on the standard electroweak parameters a,6w, Mz ■ The factor 1/Nc 
comes from the average over initial state colours. 
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We have also defined the invariants^ 



h = -2pi.h t2 = -2p2.li , (3.2) 

while is the invariant mass of the lepton pair which we fix. The component tf + t2 is 
the parity conserving piece also present in the case of the virtual photon process while the 
^1 ~ ^2 component is related to the parity violating piece of the electroweak coupling and 
hence absent for the photon case. 

We shall study the integrated cross-section which is directly related to the number of 
events below some fixed value of ay 

from which the distribution in ax can be obtained by differentiation and the dependence 
of S on will be henceforth implied. 

At the Born level the px of the lepton pair and hence the ax vanishes so that the 
full Born contribution, evaluated at fixed mass contributes to the cross-section in 
eq. ([3.31) . Evaluating this quantity is straightforward and further explanation is available 
in appendix The result we obtain is 

sW(aT) = g— 4r^ / dx, / dx2 [fq{xi)Mx2) +q^q]6{M^- sx,X2) ^ S^") , 
3vr I\c Jo Jo 

(3.4) 

where fq{xi) and fq{x2) denote parton distribution functions and xi, X2 the momentum 
fractions carried by the incoming partons. We also do not explicitly indicate a sum over 
quark/anti-quark flavours q which should be understood. In writing the above result we 
have for simplicity integrated inclusively over the lepton rapidities, which results in the 
fact that the parity violating contribution proportional to Bi Bq averages to zero. It is 
straightforward to adapt our results to include for instance acceptance cuts when final 
experimental data becomes available. 

We now derive the QCD corrections to leading order in as with the aim of identifying 
logarithmically enhanced terms in a-r to the integrated cross-section defined in eq. (p.3|). 
To this end we need to consider the process pi + P2 = h + h + k where A: is a final state 
parton emission as well as O (a^) virtual corrections to the Drell-Yan process. 

Let us focus first on the real emission contribution. We need to compute the quantity 

E«(aT)= / dxi dx2 \[fqixi)fg{x2) + q^q]t^A\aT) 

Jo Jo ^ (3.5) 

+ [{fq{xi) + fq{xi)) fg{x2) + q,q^g] sg^(aT)} , 
where the partonic quantities S^j^ which give the 0{as) contribution read 

^^Har) = j d^{hM,k)Mi{hM,k)5 (M^ - 2/1./2) e(aT - h\sin(k\) , (3.6) 



*The quantities ti and t2 were labelled as ii, £2 while h and I2 were labelled ki and k2 in ref. [23] 
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where the index i runs over the contributing subprocesses at this order, i.e. i = A/C 
denotes the annihilation (Drell-Yan)/Compton subprocesses while A4f is the appropriate 
squared matrix element the explicit form of which we include in appendix ^ We also need 
to carry out the integration over the three-body Lorentz invariant phase space since in 
addition to the final state lepton four-momenta li, I2 we also have a final state emitted 
parton k. We have introduced a delta function constraint that indicates we are working 
at fixed invariant mass of the lepton pair 2I1.I2 = M"^. Additionally in order to compute 
the integrated ax cross-section eq. (|3.3| ) , we need to restrict the additional parton emission 
k such that we are studying events below some value of a^. Recalling, from the previous 
section, that the value of this quantity generated by a gluon with transverse momentum 
kt and angle with the lepton axis (j) is kt \ sin0| we arrive at the step function in the above 
equation.^ We then fold the parton level result with parton distribution functions precisely 
as for the Born level result Yj^^^ reported above. 

After integrating over all lepton variables, accounting for virtual corrections and re- 
taining only singular terms in the limit kt ^ (which are the source of logarithms in ay), 
as detailed in appendix we arrive at the result for the annihilation contribution 

X dxi I dX2 [fq{xi)fq{x2) + q -H- q] 6 (M^ - SXlX2z) X 

Jo Jo 

x/ —Cf—^ ^ ' :e VA sin0 --V , 3.7 
Jo 27r 27r ^(1 - zf - 4zA V ' ' M/ ' ^ ' 

while that for the Compton subprocess reads 

dxi / dx2 [{fq{xi) + fq{xi)) fg{x2) +q,q^ g] 6 (M^ - SXiX2z) X 

Jo 2ir "21 y'(l - z)2 - 43A \ ' ^' Ml ^ ' 

Note that the above equations involve the step function constraint Q {kt \ sin 01 — ut) which 
represents the fact that the number of events with kt \ sin(/)| < ot is equal to the total rate 
minus the events with A;t|sin0| > ay. Since the total rate is a number independent of 
ut, we can simply compute the events with fct|sin(/)| > aj- to obtain the logarithmic ot 
dependence, which is what we have done above. 

In the above equations we have also parametrised the integral over the gluon momen- 
tum via the rescaled transverse momentum A = A;|/M^, the azimuthal angle 4> and z where 
in the collinear limit 1 — z is just the fraction of the parent partons energy carried off by 
the radiated gluon. We also have as usual Cp = 4/3, Tji = 1/2, = g'^/A-K. 



^As we stated previously this approximation is sufScient up to terms that vanish as ar — > 0, which we 
do not compute here. 
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The above results are sufficient to obtain the logarithmic structure in ut and compare 
it to the corresponding result for the Z boson pT distribution. In this respect we note that 
the only difference between the results reported immediately above and those for the pr 
case are the | sin cf)] terms in the step function constraints above. While at the leading order 
these will essentially just be a matter of detail we shall see that the sin cf) dependence has 
an important role to play in the shape of the rcsummcd spectrum. 

To complete the calculations one proceeds as in the Z boson px case and hence we take 
the moments with respect to the standard Drell-Yan variable r = thereby defining 



Jo 



cirr^-'S(ar), (3.9) 

which can be expressed as a sum over the moment space annihilation and Compton terms 
t{N, ar) = tA{N, ar) + tc {N, ar). 

The Born level Drell-Yan contribution can then be expressed in moment space as 



= Fa{N) , (3.10) 



(3.11) 



9 AiAg 

where Fa{N) denotes the moment integrals of the parton distribution functions 

Fa{N) = I dxi ,xf / dX2 x^ [/,(xi) fqix2) + Q ^ q] 

Jo Jo 
= fg{N)fg{N) + q^q, 

where we introduced f{N), the moments of the parton distributions. 
Likewise the 0{as) annihilation contribution can be expressed as 

where e = ar/M is a dimensionless version of the ot variable. 
Performing the integrals over z and A we obtain the result 



f:^')(iv,aT) = -s(°)(iv) 



e 

27r 



2Ci?as f ^d(j) 2 I sin 01 3, Isini 
H / — In In 



(3.13) 



where we introduced the quark anomalous dimension 



7<,,(7V) = Cp f dz (z^ - 1) . (3.14) 

Jo 1 — z 

Notice the proportionality of the above result to the Born level result which is a consequence 
of the collinear origin of logarithmic terms. 
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We have not integrated over the variable (j) as yet in order to make the hnk to results 
for the pt distribution. To obtain the 0{as) integrated cross-section for the px case the 
same formulae as reported above apply but one replaces | sin ^\ by unity while e would 
denote px/M. The (j) integral is then trivial and can be replaced by unity. For the ot 
variable on performing the (j) integral we use the results 



27r 



27r 



In I sin ( 



2lT 



27r 



■ In I sin ( 



ln^2 + 



In 2, 



vr 



12' 



(3.15) 
(3.16) 



to obtain 



t^l\N,aT) 



2f7.(iV)ln- 



2CFas 



TT 



In' 



2e 



5i„l 

2 2e 



+ 



Cpas vr 

2tt y 



21 



(3.17) 



The corresponding result for the QCD Compton process is purely single logarithmic 
and reads 



sg^(iV,ar) 



FciN){2^^gg{N) I ^In 



2tt 



G AiAq „ I 



sm ( 



where 



(3.18) 



(3.19) 



and Fc{N) is the moment integral of the relevant combination of parton density functions 

Fc{N) = dxi /' dx2 x^ [(/g(xi) + fg{xi)) fg{x2) +q,q^g] 
Jo Jo 



f^{N)+fg{N))fg{N) + q,q^g. 



(3.20) 



In our final results, eqs. ( pTl) and (|3T8D , we have neglected constant terms that are 
identical to those for the Drell-Yan pT distribution computed for instance in [24]. 

We note that the logarithms found here, both in the Drell-Yan and Compton contribu- 
tions, are the same as those for the px variable with the replacement e ^ 2e. In other words 
as far as the logarithmic dependence is concerned we obtain that the result for the cross- 
section for events with ax < eM is the same as the result for the variable pt/2 < eM. The 
only other effect, at this order, of the | sin0| term is to generate a constant term '^f °^ ^ 
reported above. Thus to leading order in Us we have simply 



S«(ar)| 



(PL 
2 



\pT/2=eM 



as 

27r 3 



(3.21) 



In writing the above we have returned to r space by inverting the Mellin transform so as 
to obtain the result in terms of the factor T,^^^ rather than T,^^\N). 
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The result above can be verified by using a fixed-order program such as MCFM. One 
can obtain the results for the integrated cross-sections for ut and pt/2 and the difference 
between them should be a constant with the value reported above. This is indeed the 
case, as one can see from the plot in figure ^, where the difference in eq. (3.21) generated 
using the numerical fixed-order program MCFM [21], divided by the Born cross section 
is plotted against L = ln(e). The results from MCFM agree with our expectation 
( |3.21 ) . In order to show the smoothest curve we have taken the case where the Z decay has 
been treated fully inclusively (i.e. we have not placed rapidity cuts) and a narrow width 
approximation eventually employed but we have checked our results agree with MCFM for 
arbitrary cuts on lepton rapidities. Having carried out the fixed-order computation, which 




-0.09 " ' ' ' ' ' 

-10 -8 -6 -4 -2 

L = ln(e) 



Figure 2: The difference between the integrated distributions for ar and pt/'2- Here we have used 
the CTEQ6M pdf set [25] and both factorisation and renormalisation scales and iir have been 
fixed at the Z boson mass Mz- 

serves to illustrate some important points, we shall shift our attention to the resummation 
of logarithms to all orders. 



4. Resummed Results 



Here we shall carry out the resummation of the large logarithms in the ratio of two scales M 
and ax which become disparate at small ax, ax <^ M. We already derived the dependence 
of the ot on multiple soft and/or collinear emissions in the preceding section and hence in 
order to carry out the resummation we next need to address the dynamics of multiple low 



-10- 



kf emissions. We shall first treat only the Drell-Yan process and later specify the role of 
the QCD Compton production process. 

We shall study as before the integrated cross-section representing the number of events 
below some fixed value of a-r, defined in eq. ( |3.3D , from which one can obtain the ax distri- 
bution by differentiating with respect to a-j-. Also as we emphasised in the previous section 
we are working at fixed invariant mass of the lepton pair purely as an illustrative example 
and we can straightforwardly adapt our calculations to take into account experimental cuts 
on for instance lepton rapidities, which in any case do not affect the resummation. 

We consider again the incoming partons as carrying momentum fractions xi and X2 
of the incoming hadrons which means that at Born level where they annihilate to form 
the lepton pair via virtual Z production we have simply = s = sxiX2, where the 
Mandelstam invariant s denotes the partonic centre of mass energy. Beyond the Born 
level one has to take account of gluon radiation and to this end we introduce as in the 
previous section the quantity z = M'^/s such that 1 — z represents the fractional energy 
loss of the incoming partons due to the radiation of collinear gluons. Thus in the limit 
2; — )• 1 one is probing soft and collinear radiation while away from 2; = 1 we will be dealing 
with energetic collinear emission. We note here that for the purpose of generating the 
logarithms we resum we do not have to examine large-angle radiation and the collinear 
limit is sufficient as for the usual pT distribution. In fact since the ay resummation we aim 
to carry out shares several common features with the well-known pT distribution we shall 
only sketch the resummation concentrating instead on features of the a-r which lead to 
differences from the pT variable. For a recent detailed justification of the approximations 
that lead to NLL resummation for the px case as well as for other variables the reader is 
referred to ref. [26]. 

We work in the centre-of-mass frame of the colliding partons and in moment space 
where we take moments with respect to r = M^/s of the cross-section in eq. ( |3.3| ) as in 
the fixed-order calculations we carried out. Taking moments enables us to write for the 
emission of multiple collinear and optionally soft gluons 

t{N,aT) = t''^\N)WN{aT), (4.1) 



where 'E^'^\N) is the Born level result in eq. ( |3.10 ). The effects of multiple collinear (and 



optionally soft) gluon emission from the incoming projectiles are included in the function 
Wn which can be expressed to next-to-leading logarithmic (NLL) in the standard factorised 
form 



00 ^ n „ 



n=0 1=1 



x.f2C^^(i±^)e aT-|E-(^^)lh (4-2) 



where 1 — Zi denotes the fraction of momentum carried away by emission of a quasi-collinear 
gluon i from the incoming hard projectile so that M"^ /s = z = W - Zi and ku is the transverse 
momentum of gluon i with respect to the hard emitting incoming partons. 
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In writing the above results we have used an independent emission approximation, 
vahd to NLL accuracy, where the emission probabihty for n cohinear gluons is merely 
the product of single gluon emission probabilities, which factorise from the Born level 
production of the hard lepton pair.^ The single gluon emission probability to the same 
NLL accuracy is given by the leading order splitting function for the splitting of a quark 
to a quasi-coUinear quark and gluon (weighted by the running strong coupling), 

n,(.)^^C.=^i±f, ,4.3) 

with Os defined in the CMW scheme [29] . We have inserted a factor of two to take account 
of the fact that there are two hard incoming partons which independently emit collinear 
gluons. We have also taken care of the constraint on real gluon emission, imposed by the 
requirement that the sum of the components of the ku normal to the axis in eq. ( ^.1| ) 
(denoted by v{ki) = ku siiKpi) is less than ax- We have integrated over the leptons, holding 
the invariant mass M fixed, and taken moments to obtain the full zeroth order Drell-Yan 
result T,^^\N), which multiplies the function Wn containing all-order radiative effects. 

All of the above arguments would also apply to the case of the pT variable. Thus 
while the dynamics of multiple soft/collinear emission is treated exactly as for the pt 
resummation the difference between the pT and our resummation arises purely due to the 
different form of the argument of the step function restricting multiple real emission. Thus 
while for the pT variable the phase space constraint involves a two-dimensional vector sum 
B (^pt — I Yli ^ti\j , in the present case we have a one dimensional sum of the components 
of the gluon kt normal to the lepton thrust axis, v{ki) = kusiiKpi. One encounters such a 
one dimensional sum also in cases such as azimuthal correlations in DIS [20, 30] and the 
resummation of the pt difference between jets in dijet production [22]. It is this difference 
that will be responsible for different features of the ay distribution as we shall further 
clarify below. The relationship between azimuthal correlations and the a-r is no surprise 
since the ax variable is proportional to vr — Acpu, the deviation of the azimuthal interval 
between the leptons from its Born value vr. 

In order to further simplify eq. ( [4. 2]) we also factorise the phase space constraint using 
a Fourier representation of the step function [30] 

@(^aT-\Y,v{h)\^=^ 1^ jsm{baT)lle"^<^''^\ (4.4) 

Note the presence of the sin(6aT) function which is a consequence of addressing a one 
dimensional sum as opposed to the Bessel function Ji one encounters in resummation 
of the Pt cross-section. With both the multiple emission probability and phase space 
factorised as above it is easy to carry out the infinite sum in eq. ([4. 2)) which yields 

WT\aT) = - sHbar) e^^'^'^ , (4.5) 



^This approximation is invalid for situations when one is examining soft radiation in a limited angular 
interval away from hard emitting particles [27, 28] , which is not the case here. 
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with the exponentiated real gluon emission contribution 

^reai(.) ^ [dz^^z^2 ( \±^] e^^^(^) G f 1 - . - ^ V (4.6) 

J 2tt 2Tr \1-z J \ MJ ^ ' 

The kinematic hmit on the z integration is set in such a way that one correctly accounts 
for soft large angle emissions. 

Next we include all-order virtual corrections which straightforwardly exponentiate in 
the soft-coUinear limit to yield finally 

WN{aT) = - ^siu{haT)e~^^^\ (4.7) 

where 

^ ' J kf 277 2Tr \1- z J 



(.^e^H^)-l)e(l-.-^), (4.8) 



X \^z e 

where it should be clear that the term corresponding to the —1 added to the real contri- 
bution e''^'"^^^ corresponds to the virtual corrections. Note that the virtual corrections 
are naturally independent of both Fourier and Mellin variables h and N respectively since 
they do not change the longitudinal or transverse momentum of the incoming partons 
and hence exponentiate directly. We are thus left to analyse R{h) the "radiator" up to 
single-logarithmic accuracy. 

4.1 The resummed exponent 

Here we shall evaluate the function Rib) representing the resummed exponent to the re- 
quired accuracy. We shall first explicitly introduce a factorisation scale render the 
integrals over kt finite. Later we will be able to take the Qo — ^ limit. Thus one consid- 
ers all emissions with transverse momenta below Qq to be included in the pdfs which are 
defined at scale Qq such that the factor Ti^'^\N) reads 

^^'\N) = ^^Fj,{N,Ql), (4.9) 

with 

Fa{N, Ql) = C dxixf C dx2X^ [/,(xi, Ql) + q^ q] 

Jo Jo (4.10) 

Thus while in the fixed-order calculation of the previous section the pdfs could be treated 
as bare scale independent quantities for the resummed calculation we start with full pdfs 

evaluated at an arbitrary (perturbativc) factorisation scale. The kt integration in the 
perturbative radiator should now be performed with scale kt > Qq. We next follow the 
method of ref. [31] where an essentially identical integral was performed for the radiator. 
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First, following the method of rcf. [31] we change the argument of the pdfs from Qq to 
the correct hard scale of the problem, the pair invariant mass M, via DGLAP evolution. 
To be precise we use for the quark distribution 

UN,Ql) = f,{N,M')e-^'''o -4 - '""^ \ (4.11) 

and likewise for the anti-quark distribution where ^qq{N) is the standard quark anomalous 
dimension matrix. Note that we have not yet considered the QCD Compton scattering 
process and the corresponding evolution of the quark pdf from incoming gluons via the jqg 
anomalous dimension matrix, which we shall include in the final result by using the full 
pdf evolution rather than the simplified form reported immediately above. Carrying out 
the above step results in a modified radiator such that one now has 

X 



X "rd.i±i!fi-.^e^M^))efi-.-^V^ 

Jo l-z \ J \ Mj Cf 



(4.12) 



Using the definition of the anomalous dimension ^qq{N) we can write the above (see for 
instance ref. [31]) as 

^ ^ ^7^2 27, Jo 21, Jo l-z\ ) \ MJ ' 

(4.13) 

where in arriving at the last equation wc neglected terms of O (kt/M). 
To NLL accuracy we can further make the approximation [30] 

1 _ e'^^W ^ (^kt\ sm(j)\ - , (4.14) 

where we used the fact that in the present case v{ki) = kt sin (f)i with b = be^'^. Thus one 
gets for the radiator R{b) = R{b) with 

M2 ,,2 tu2 



xe(^fci|sin0|-i)e(l-^-^) . (4.15) 

We now evaluate the radiator to the required accuracy. Performing the z integration to 
NLL accuracy, and neglecting terms of relative order kt/M, one arrives at 

= f""'4^s(k,lsi.^l-l]LMl-U'-^) . (4.16) 



2tt Jo k1 27r V bj \ kf 2 Cp 

where in the last line we took Qq^O since the kt integral is now cut-off by the step function 
[fetjsin^l — 1/6] and chose to perform the cp integration at the end. We note that the 
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only difference between the radiator above and the standard distribution is the factor 
I sin(/>| multiplying fcj in the step function condition above. This will result in an additional 
single- logarithmic contribution not present in the px resummation results. 

In order to deal with the (j) dependence to NLL accuracy we expand (as in ref. [22]) 
eq. ( P^ about I sin (f)\ = 1 m. powers of In | sin (f)\ to obtain 



R{b) = 2Cf 



M2 



27r 



In 



M2 



3 lUN) 
2 Cf 



+ 



dR(b) 



2n 



d<P . . . 
— In sm ( 
27r ' 



^, , (4.17) 

dln{bM) 70 2^ 

where we have neglected higher derivatives of R as they will contribute only beyond NLL 
accuracy. Moreover in evaluating dR{b)/d\n{bM) we can replace R by its leading loga- 
rithmic form i?LL(^) since logarithmic derivatives of any next-to- leading logarithmic pieces 
of R{b) will give only NNLL terms that are beyond our accuracy. The first term on the 
RHS of the above equation is in fact just the radiator we would get for resummation of 
the Z boson px distribution which contains both leading and next-to-leading logarithmic 
terms. The second term on the RHS accounts for the dependence of the problem and is 
purely next-to-leading logarithmic in nature since it contains the logarithmic derivative of 
i?LL(^)- Thus we need to evaluate the first term on the RHS of eq. ( [4.16|) and then isolate 
its leading- logarithmic piece to compute the second term on the RHS. 

We carry out the integral over kt to NLL accuracy using standard techniques [32], 
i.e. we change the coupling from the CMW to the MS scheme 



a 



CMW 



a 



MS 



(kl 



1 + K 



a 



MS(fc2) 



2tt 



A- = C.I 5!-^ 

'18 6 



and use a two-loop running coupling 



as{kl) 



l-p 



l-a.(M^)fMl^ 



a,(M2)/3oln 



k? 



where the beta function coefficients are defined as 



/3o 



IICa - 2nf 



/3i 



17C\-bCAnf -^Cpuf 



(4.18) 



(4.19) 



(4.20) 



127r ' 247r2 
and as(M^) is a shorthand for a^^(M^). We then obtain the usual expression 

R{b) = LgiiasL) + (72(0^^) , (4.21) 

with L = In (JPM^) . The functions gi and §2 are then the leading and next-to-leading 
logarithmic functions which have the following detailed form 

Cf 



51(A) 
52(A) 



3Cf 
''2^ 



[-A-ln(l-A)] , 

ln(l-A)-^ln(l-A) + ^A(-ln2) 



7r/3o 



KCj.[A + (l-A)ln(l- A)] _ CfPi 
27r2/32(i _ A) ^ 



7r/3o 1 - A 
A + ln(l-A) 1 

1 - A 2 ^ 



A) 



(4.22) 



(4.23) 
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with A = as{M^) PqL. 

Let us comment on the origin of various terms. The leading logarithmic function 51(A) 
arises from soft and collinear emission integrated over the phase space with a running 
coupling (to be precise the one-loop running of the coupling is sufficient to give us gi). It 
is identical to the corresponding function for Z boson pT resummation and at this level 
the ot and pT variables do not differ. The function §2 embodies hard-collinear radiation 
(and hence the appearance of the quark anomalous dimension ^qq{N)) as well as the two- 
loop running of the coupling and the change to the MS scheme from the CMW scheme 
which gives rise to the piece proportional to K. It is also the same as the corresponding 
function for Z boson pT resummation except for the additional single-logarithmic term 
^^Y3j(— ln2) which arises from the (j) dependence of the problem. In other words one 
has explicitly 



d\a{hM) 27r ' 51n(6M)' ' 7r/3o 1 - A 

However, this term, within NLL accuracy, can be absorbed in the radiator with a change 
in the definition of its argument h: 

and R(b/2) = R{be^^ /2) is precisely the radiator for the Z boson px distribution (see 
e.g. ref. [19]). As a final step we can use the anomalous dimension matrix jqq{N) and 
the corresponding contribution "yqg{N) from Compton scattering which we have for brevity 
avoided treating, to evolve the pdfs from scale to scale (2/6)^ precisely as in the case 
of the pt variable. After absorption of the dependent piece of the radiator into a change 
of scale of the pdfs it is trivial to invert the Mellin transform to go from A^ space to r 
space. We can thus schematically write the result for the ax cross-section defined in 
resummed to NLL accuracy as 

,APAiAq2 db 



, ^M^AiAq2 db . , , 
S AT = G- / -r smibar) exp -5 6/2 x 

StT Nc 7T Jo b 



Sm =2C,I _ In ^ - ) . (4.27) 



dxi / dx2 [fq (xi, (2/6)2) (x2, {2/bf) +q^q] 6{M^ - xiX2s) . (4.26) 
Jo 

The function 5(6/2) is just the radiator i?(6) without the A^ dependent anomalous dimen- 
sion terms which have been absorbed into the pdfs: 

dfcf a,(fc|) / 3^ 
/(2/6)2 fcf 27r \ 2^ 

As we just pointed out this function coincides with the corresponding function in the 
Pt case. All differences between ax and pT thus arise from the fact that one has to 
convolute the resummed exponent and pdfs with a sine function representing the constraint 
on a single component of kt, rather than a Bessel function representing a constraint on 
both components of the kf. Having achieved the resummation we shall next expand our 
resummed result to order and compare its expansion to fixed-order results to non- 
trivially test the resummation. 



-16- 



5. Comparison to fixed-order results 



In order to non-trivially test the resummation we have the option of expanding the results 
to order (i.e. up to two-loop corrections to the Born level) and testing the logarithmic 
structure against that emerging from fixed-order calculations. Since the results for the pT 
distribution are already well-known and since many terms are common to the ut and pT 
resummed results it is most economical to provide a prediction for the difference between 
the ct and pT variables. To be precise we already identified a leading-order result for the 
difference between cross-sections involving ut and in eq. ( p. 21 ). In this section we 



shall derive this difference at NLO level and compare to fixed-order estimates. 

Let us consider the resummed results for the ay and pt/2 cases. We remind the reader 
of the well-known result for the pT variable by expressing the integrated cross-section for 
events with below a fixed value eM: 

S (pt/2) \p^/2=eM = (l + Ci(iV) ^) s(°)(iV) J^^db2Me Ji{h2Me) e"''(^) , (5.1) 

The above result is expressed in moment space and we have additionally provided a mul- 
tiplicative coefficient function (1 -|- Ci(A^) as/(27r)), so that the result accounts also for 
constant terms at leading order. This form of the resummation is correct up to NNLL 
accuracy in the cross-section whereas the pure resummed result without the multiplicative 
constant piece is correct to NLL accuracy in the resummed exponent [32]. Thus with the 
constant Ci in place the resummation should guarantee at order a1 terms varying as a^L^, 
a'lLp' as well as the a^L^ term which partially originates from a "cross-talk" between the 
agCi term and the a^L^ term in the expansion of the exponent.^ 
The equivalent result for the ax variable is 

S(aT)|a,=.M= (l + C'i(iV)^)sW(iV)^^"* sin(6Me)e-^(^) , (5.2) 

where Ci is the constant for the ay variable. We shall first expand the resummation to 
order as and consider the difference between ax and pt/2. First we express the radiator 
in the standard notation [32] 

oo n+1 

-R (6/2) = ^ ^ a/L- , L = In {b'M^/A) , (5.3) 

n=0 m=0 

with ds = as/27r. Having done so we expand the resummed exponent so that to order 
we can write for the pT variable 



poo 

S {pt/2) |p^/2=,M = (1 + Ci(iV) ds) S(0)(iV) / db2Me Ji{b2Me) x 

JO 

X (1 + GiidsL + Gi2dsL'^ + 0{al)) , (5.4) 

^This form of the result we use is an oversimplification since we consider only the piece of the 0{aa) 
constant which is associated to the annihilation channel. In principle one should also include the constant 
arising from the Compton channel but this is identical to the corresponding constant for the pr variable [24] 
and it is straightforward to show that its effects cancel to the accuracy we need for the result we derive 
below for the difference between ar and pt/2 variables. 
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where we replaced the resummed exponent by its expansion to order as- 
Carrying out the b integral yields 



S {pt/2) 



PT/2=eM 



(1 + Ci{N)as) s(°)(iV) ( 1 + Giia.ln 



" 1 " 




" 1 ■ 




4e2 


+ Gi2as In^ 


4e2 


) 



(5.5) 



where for the moment we do not insert the explicit forms of the Gnm coefficients. 
Repeating the exercise for the ay variable one obtains 



^{aT)\ar=eM = (l + Ci{N)ds) SW(iV)x 

1 + GiiOs In 



1 ■ 




" 1 " 


4g2 


+ Guas In^ 


4e2_ 



+ Gi2 Os 



IT 



(5.6) 



where we labelled the constant piece as Ci to distinguish it from that for the pT variable. 
Constructing the difference at 0{as) between the ut and pt/2 variables we find that all 
the logarithms cancel and we obtain 



^{aT)\aT=eM - S (Pt/2) \pj,/2=eM 



a. 



IT 



Ci{N)-Ci{N) + Gi2Y 



(5.7) 



The value of the resummation coefficient G12 can be obtained from eq. ( 4.23 ) by 
expanding the result in powers of A from which we find G12 = —Gp- Comparing this result 
with our explicit leading order calculation eq. ( |3.21| ) we find that Gi = Ci. 

Next we carry out the expansion of our resummation to order and construct the 
difference from pt/2 at this order. We shall then compare our expectation with MCFM. 
Expanding the radiator to order one gets 



-R 



+ a 



1 + ds {GuL + GuL'') 
Gl,L^ 



+ G22^^ + GiiGi2L^^ 



G23L^ + 



(5.8) 



Retaining only the 0{a'^) terms we can once again carry out the b space integrals as 
before and in particular the new integrals that appear at this order are 



Tip) 

ax 



Tip) 
PT 



db 



sm{bMe) In^ b 



db2MeJi{b2M()W b' 



,M2 



(5.9) 



with p = 3,4. Carrying out the above integrals with p = 3 is straightforward and the 
difference between the integrals for the pT and ay case with p = 3 produces only an 
a^lnl/e^ term apart from constant pieces. Such terms are beyond the accuracy of our 
resummation which ought to guarantee only terms as singular as In^ 1 /e^ in the cross- 
section and hence to our accuracy there will be no contribution for p = 3 for the difference 
between ar and pt/2. 
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Figure 3: The difference between D{aT) and D{pt/2), defined in eq. ( 5.1f ) witfi tlie subtraction 
of tlie computed logarithmic enhanced term in eq. (|5.10|) . 



The situation changes when we consider the p = 4 integrals. All relevant logarithms 
cancel between ax and pt/2 once again except a term varying as aj In^ ^. To be precise 
considering only the order a.g terms one obtains 



j:{aT)\aT=eM " S (pt/2) 



\pt/2=€M 



^ X 27r2a2 ln2 1 x sW(7V) + Oia^L) 



vr^ C| In^ 



s(°)(iv). 



(5.10) 



To convert the result above back into r space from Mellin space is straightforward as 
one just inverts the Mellin transform for T,^^^ to yield the Born- level quantity T,^^^ as the 
multiplicative factor. 

Once again the above result can be tested against the results from MCFM. We con- 
sider the difference in the differential distributions (derivative with respect to Ine of the 
appropriate integrated cross-sections) for ax and as a function of In e 



Diar) 



-D{pt/2) 



1 



\pT/2=eM 



din e 



\aT=eM 



dine 



\pT/2=eM 



(5.11) 



Our prediction for this difference can be obtained by taking the derivative with respect to 
Ine of the RHS of eq. ( ^.10 ). Subtracting this prediction from the MCFM results should 
yield at most constant terms arising from the logarithmic derivative of formally subleading 
a1 Ine terms. That this is the case can be seen from fig. ^ where we note that at sufficiently 
small values of e the difference between MCFM and our prediction tends to a constant. 
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6. Discussion and conclusions 



Before concluding we should comment on the resummed result eq. ( [4.26| ). First we note 
that there is the usual issue that is involved with b space resummation of the large and 
small b behaviour of the integrand in that the resummed exponent diverges in both limits. 
The small b region is conjugate to the large kt regime which is beyond the jurisdiction of 
our resummation. At sufficiently large b on the other hand we run into non-perturbative 
effects to do with the Landau pole in the running coupling. These issues can be resolved 
by modifying the radiator such that the perturbative resummation is not impacted. For 
instance the strategy adopted in ref. [33] was to replace the resummation variable b by 
another variable b* which coincides with b in the large b limit but at small b ensures that 
the radiator goes smoothly to zero. Likewise to regulate the Landau pole a cut-off was 
placed in the large b region of integration in the vicinity of the Landau pole and it was 
checked that varying the position of the cut-off had no impact on the resummation. Other 
prescriptions can be found for instance in [34]. 

As far as the behaviour of the resummed cross-section and consequently the corre- 
sponding differential distribution is concerned the difference from the pT distribution is 
solely due to the convolution of the resummed b space function with the sin(6) function 
rather than a Bessel function. As was explained in detail in ref. [30] the result of convo- 
lution with a sine function produces a distribution that does not have a Sudakov peak. 
The physical reason for this is that a small value of qt can be obtained by two compet- 
ing mechanisms. One mechanism is Sudakov suppression of gluon radiation and this is 
encapsulated to NLL accuracy by the resummed exponent. The other mechanism is the 
vectorial cancellation of contributions from arbitrarily hard emissions which in this case 
involves cancellation only of a single component of kt transverse to the lepton axis. This 
mechanism is represented by the presence of the sine function while a two-dimensional 
constraint such as that for the pT variable is represented by a Bessel function. In the case 
of one dimensional cancellation such as for the ax as well as for instance for the dijet A<j) 
variable [20] the cancellation mechanism dominates the Sudakov suppression mechanism 
before the formation of the Sudakov peak while for the pT case the vectorial cancellation 
sets in as the dominant mechanism after the formation of the Sudakov peak. Thus for 
the qt distribution one sees no Sudakov peak but the distribution rises monotonically to 
a constant value as predicted by eq. ( [4.26 ). 

To conclude, in this paper we have carried out a theoretical study based of a variable, ot 
proposed in ref. [1] as an accurate probe of the low pT region of the Z boson pT distribution. 
Having accurate data on the ot well into the low ut domain will be invaluable in pinning 
down models of the non-perturbative intrinsic kt and may lead to firmer conclusions on 
aspects such as small-x broadening of pr distributions [8] than have been reached at present 
with the pt variable. In this respect it may also be of interest to examine theoretically 
the power corrections to the aj- distribution along the same lines as for the pT case [35] 
and hence to examine theoretically whether the ax and pT ought to have identical non- 
perturbative behaviour. This is once again work in progress. 

Before any such conclusions can be arrived at however, it is of vital importance to have 
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as accurate a perturbative prediction as possible to avoid misattributing missing perturba- 
tive effects to otfier sources. Tlie most accurate perturbative prediction one can envisage 
for the ax case is one where resummation of large logarithms in ax is supplemented by 
matching to fixed-order corrections up to the two-loop level. In this paper we have carried 
out the first step by resumming to NLL accuracy the ax distribution and checking our 
resummation by comparing to the logarithmic structure in exact fixed order calculations. 
We envisage that it should be possible to actually extend the accuracy of the resummed 
calculation to the NNLL level which has already been achieved for the px variable [19] and 
this is an avenue for future development. In any case our current prediction matched to 
fixed-order estimates from MCFM should already enable accurate phenomenological inves- 
tigation alongside forthcoming Tevatron data [1]. As part of an article in progress we plan 
to carry out the matching and a detailed phenomenological study of the ax distribution 
which we anticipate will shed more light on issues relevant to physics at the LHC in the 
near future. 
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A. Born level result 

Here we explicitly compute the leading-order and the real part of the next-to-leading order 
contribution to the ax integrated cross section defined in eq. (^^). 
At leading order Ti^^\ax) is just the Born cross section: 

S(0)(a5n) = I dxi [ dX2 [fq{xi)fq{x2) + q ^ q] X 
Jo Jo 

X J d^{liM)MlY{hM)5{M^ -"^h-h) , (A.l) 

where A^^Y Born matrix element reported in eq. ( |3.lD . To this end we look at the 

Lorentz-invariant phase-space which can be written as 

where we included in addition to the usual two-body phase space a delta function cor- 
responding to holding the invariant mass of the lepton-pair at and s is the partonic 
centre of mass energy squared s = s xiX2- 
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We parameterise the four vectors of the incoming partons and outgoing leptons as 
below (in the lab frame) 



pi = ^xi(l,0,0,l), (A.3) 

P2 = ^X2 (1,0, 0,-1), 
^1 = (cosh y, 1, 0, sinh y) , 

with I2 being fixed by the momentum conserving delta function. 

In the above ^/s denotes the centre of mass energy of the incoming hadrons while 
xi and X2 are momentum fractions carried by partons pi and p2 of the parent hadron 
momenta, while It and y are the transverse momentum and rapidity of the lepton with 
respect to the beam axis and we work in the limit of vanishing lepton masses. In these 
terms we can express eq. ( |A.2| ) as (after integrating over I2 using the momentum conserving 
delta function) 

ldHh,l2) = ^J'-^S{{p,+p2-hf), (A.4) 

where we have carried out an irrelevant integration over lepton azimuth. Note that the 
factor 6 [{pi + P2 — h)'^) arises from the vanishing invariant mass of lepton 12- 

In order to obtain the full Born result we need to fold the above phase space with the 
parton distribution functions^ and the squared matrix element for the Drell-Yan process 
to finally obtain 



sW(M2)= f\xif{xi) C 
Jo Jo 



X 

2s 



dX2 /(X2) X 

1 f Ix dlx dy 



j ^^4^'^^ '^ ^1^2 +tl+ t2) 6 (M2 - XIX2S) M^y. (A.5) 



where we used (pi + P2 — lif' = sxiX2 + ii + ^2- 

We next evaluate the squared matrix element A^^y ( |3-1| ) in terms of the phase 

space integration variables: 

ti = -2pi.li = —^/sxllTe~^ , t2 = -\/sx2lTe^ ■ (A-6) 



Inserting these values of ti and t2 in eq. (|A.5D we use the constraint 

5{sxiX2 + ti+t2) = 5 {sxiX2 - y/slr {x2 + xi e~^)) , (A. 7) 

to carry out the integration over It which gives 

C dxi f{xi) C dX2 f{x2) 5 (M2 - XIX2S) ^ ^ Mly , (A.8) 

8vrs Jo Jo {x2ey + xie-y) 



*In order to avoid excessive notation we do not explicitly indicate the sum over incoming parton flavours 
which should be understood. 
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tl = xle-^y- ^ — , (A.9) 



where in evaluating Ad^y '^^^ needs to use It = \/s xiX2/{x2e^ + xie ^) , which yields 
using ( |A.6D 

{x2ey + xie~y) 
M 

{x2ey + xie^y 
Using the above to evaluate A^^y ( p-l] ) we obtain 



tl — x\e^y- — ——2 . 



s(o)(m2) = ^ / dxifixi) dx2f{x2)6{M^ -xiX2s) dyF{xi,X2,y), (A.IO) 

Nc TTS Jo Jo ^ J 



where we introduced 



T (xi , X2, y) = A A + Bi B, \ ' . (A.ll) 

[X2ey + xie y) [X2ey + xie 



Integrating the angular function over rapidity over the full rapidity range^ one finds as 
expected that the parity violating component proportional to Bi Bg vanishes and the result 
is Ai Aq/{3xiX2)- Thus the final result is (using xiX2 = M^/s) 

s(0)(M2) = g— ^ / dxifixi) dx2f{x2)6{M^-XiX2s). (A.12) 

-iTT Jyc Jo Jo 

We now compute the LO QCD result by considering the emission of a gluon in the 
Drell-Yan (QCD annihilation) process as well as the contribution of the quark-gluon (QCD 
Compton) scattering process. Thus we consider the reaction pi + p2 = h + h + k where 
k is the emitted gluon in the Drell-Yan process and a quark/anti-quark for the Compton 
process. 

The squared matrix elements at this order for the annihilation and Compton processes 
are reported in appendix ^ and we shall use those results in what follows below. 
We explicitly parameterise the momentum k as below 

k = kt (cosh yfc, cos 0, sin (p, sinh y/c) . (A. 13) 



The parameterisation of the other particles four-momenta is as in the Born case eq. 
One now has to integrate the squared matrix elements over a three body final state and to 
this end we introduce the usual Mandelstam invariants 

u = -2pi.k = -y/sxihe'y'' , i= -2p2.k = -y/sx2ktey'' . (A. 14) 

The Lorentz invariant phase space is now 

/I f d^l d^l d^k 

= Ts J 2i^,2i2^„2i2^J'^'''''^ + « -k-k-k). (A.15) 

®We can straightforwardly adapt the calculation to include the experimental acceptance cuts when 
available. 
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Following the same procedure as in the Born case we perform the trivial integration 
over I2 and obtain the leading order QCD correction to the Born result ( [A.5| ) (for the 
moment we are considering just real emission terms indicated below by the label r) 



S«(M2)= V / dxifixi) [ dx2f{x2)x 

1 f Ixdlxdy d'^k 



^sUl^-.ms{l-^!i-'^J^\\Ml (A.16) 

Here the factor ^1 — _ accounts for the energy-momentum carried off by the ra- 

diated parton k while the index i = A pertains to the QCD annihilation process while i = C 
indicates the QCD Compton process. Noting that one has as before ti = —^/s xilTe~y ,t2 = 
—y/sx2lTe^ and additionally 2li.k = 2 It h {cosh{y — y^) — cos (p), we can use the con- 
straint 6 (M^ + ^1 + ^2 + 2/1. /c) to integrate over It and the value of It (and hence ti,t2) 
is thus fixed in terms of other parameters: 

It 



^/s l^xie y + X2ey — 2-^ (cosh(y — y^) — cos ( 

xie~yM^ 

xie~y + X2ey — 2-h= (cosh(y — y^) — cos (p) 



(A.17) 



with the expression for t2 the same as that for ti except that xi in the numerator of 
the above expression for ti is to be replaced by X2 e^. After integrating away the It one 
gets 



S«(M2)= V f dxif{xi) C dX2f{x2)^ 



d^k 
^2{2TTfkQ 



X ^Mj 5{m'^ - zxiX2s) . (A.18) 

s {xie-y + X2ey - ^ (cosh(y - y^) - cos 0)^ 

where we introduced z = 1 — 2{pi.k) / s — 2{p2.k) / s . 

We are now ready to integrate over the parton and lepton phase space variables. Since 
we are interested in the specific cross-section in eq. ( |3.3| ), we need to integrate over the 
phase space such that the value of the qt is below some fixed value. Further we are 
interested in the small qt logarithmic terms so that we consider the region ut/M <C 1. 

To avoid having to explicitly invoke virtual corrections we shall calculate the cross- 
section for all events above qt and subtract this from the total O (as) result T,^^\M'^) 
which can be taken from the literature [36]: 

sW(ar,M2) = ^^(M^) - ^^(aT.M^) , (A.19) 
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where we shall calculate T!i'\aT,M'^) = J^^ HiriM^'^^'T- 

Moreover since we are interested in just the soft and/or collinear logarithmic behaviour 
we can use the form of the ax in the soft/collinear limit derived in the previous section. 
Thus we evaluate the integrals in eq. (|A.1S| ) with the constraint {kt \ sva.(j)\ — ax)- In order 
to carry out the integration let us express the parton phase space in terms of rapidity y^, 
kt and (j). Thus we have 

d'^k f ktdktdykdcp 



2(27r)3A;o J 2(27r)3 

M'^\ f d(h f, r dzdA 



where we used 



(A.20) 



1 - -^e-y>' - -^eV' , (A.21) 



which follows from the definition of z and where we also introduced the dimensionless 
variable A = kl/M"^. A fixed value of z corresponds to two values of the emitted gluon 
rapidity 



y± = In 



'I- z)± {I - zY - 4:Z/S)j 



(A.22) 



2kt 

The requirement that the argument of the square root in eq. ( A.22| ) be positive sets an 



upper bound for the z integration, specifically z < 1 — 2{^J A + A^ — A). In the limit 
A — )• 0, up to corrections of order A, this bound reduces to z < 1 — 2^fK. 

Having obtained the phase space in terms of convenient variables we need to write 
the squared matrix elements in terms of the same. We first analyse the QCD annihilation 
correction and next the Compton piece. In the annihilation contribution one has singu- 
larities due to the vanishing of the invariants t and u with the \/{tu) piece contributing 
up to double logarithms due to soft and collinear radiation by either incoming parton 
and the \/t and 1/n singularities generating single logarithms. The double logarithms 
arise from low energy and large rapidity emissions (soft and collinear emissions) while the 
single- logarithms from energetic collinear emissions, hence it is the small kt limit of the 
squared matrix elements that generates the relevant logarithmic behaviour. Thus we write 
the squared matrix element M.\ in eq. ( p.lD in terms of the variables A and z and then 
find the leading small A behaviour. Specifically in the A — t- limit, considering only the 
\/i singular piece, the factor appearing in ( |A.18 ) has the following behaviour (keeping for 



now only the AiAq piece of the matrix element): 



.Ml 



(^xie~y + X2ey - ^ (cosh(y - y^) - cos ( 



16,^g^^(l + .^)^^ii^:!^±^. (A.23) 
^ Nc A ^ ^ ' s {xie-y + X2zey)^ ^ ' 
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Performing the integral over all rapidities y of the lepton, the above factor produces 
16c/2g^^ld^. The corresponding 5/ piece of the squared matrix element vanishes 
upon integration over all rapidities. 

Since the 1/u singular term, after integration over all lepton rapidities, gives us the 
same result as that arising from eq. ( A.23| ), we can write for the annihilation process (using 
eqs. ( [AISD , (|A^)) and g"^ = Airas 



.(1) 



.AlAg 



dxidx2f{xi)f{x2)^T— dz 6 (M"^ - XIX2Z s) 

^'^ Jo ^ Jo 



d±^ as 2{l + z') 
27r ^27T ^{i- zf -4zA 



e(^/A|sin,/.|-^) . (A.24) 



Following the same procedure for the Compton process one finds instead 



s ^ 



s (^xie~y + X2ey — ^ (cosh(y — y^) — cos <; 

^ .AiA^Tr w2^n 2 xle-^y + xl^^ 

^Wg g — {l-z)[z + (1 - z 

Nc A s {xi e-y + X2 z ey) 



, (A.25) 



which after integration over the rapidity y reduces to 16g^^^^^^(l — z) ^ . Thus 

we have for this piece 



I-2VA 



2lT 



M2 dA 

as {l + z) [z^ + {l-z^)] 



'•^^^ I dxidx2f{xi)f{x2)"-^ I I dz 6 {M"^ - X1X2Z s) 



27r 27r 



AzA 



G(\/A|sin,/.|-^) . (A.26) 



B. Leading Order Cross Section in the small A limit 

The matrix element squared for the QCD annihilation process from ref. [23] is (in four 
dimensions) 



M\{liM,k) = -IGg'^g^M'^x 



Nc 

^ s-2ti-M^ tl + 4 + s(ti + t2 + M^)\ , . ' 
1 + k 1 + (u o t, ti o ta) 



t 



tu 



-BlBq 



[s + 2t2 + M2) M2 (ti - t2) 



+ 



tu 



(u o ■(-)• t 



, (B.l) 
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while for the QCD Compton process, if p2 represents an incoming gluon, one has 



i-2 (h + M2) 



+BiB„ 



^ s + 2{ti+t2) ^ 2 
i si 

t ^ S + 2(ti+t2) 



2M2 (2ti +t2 + M2) 

St 



(B.2) 



where we corrected small errors (after an independent recomputation of the above) of an 
apparent typographical nature in the BiBg piece of the annihilation result. 

Here we report also the electroweak coefficient constants G, Ai, Ag, Bi, Bq for the case 
of Z boson exchange: 

g{a,ew,M\Ml) = ,i^4 ll "„^4 (m^ _ + (TzMz? ' (B.3) 

Af = a)+b), Bf = 2afbf, {f = l,q). 

All these quantities have been taken from ref. [23], where the reader can find analogous 
expressions for the case in which a virtual photon is exchanged as well. Following the 
conventions of ref. [23] , we also have 



ai = - J + sin^6lH^, 

1 2 2 

au,c = - - - sm Ow , 

1 1-2/, 

ad,s,b = ~4 + 3 ^^"^ ' 



bu,c — 



; 

1 

4 ■ 



(B.4) 



The kinematical variables u and i are defined in eq. ( [A. 14 ). In the limit of small A both 
matrix elements become collinear singular. In the annihilation subprocess this occurs when 
the emitted gluon k is collinear to either pi (corresponding to n — t- 0) or p2 (t — t- 0). The 
singularity for u — )■ occurs at positive gluon rapidity y^, correspondingly the one for t — t- 
occurs at negative y^. The matrix element for the Compton process shows only a collinear 
divergence when an outgoing quark is collinear to the incoming gluon, corresponding to 
t — )• 0. In the following we compute the approximated expression of M\ and Ai"^ in the 
collinear limit t — )• 0. The remaining collinear limit u — )• of M\ gives an identical result 
after integration over the lepton rapidity. 

Neglecting terms of relative order kt one has 

M2 



It 



/s (xie y + zx2ey) 



t 



1 



t2 



1-z' 
X2eyM^ 



(B.5) 



{xie y + zx2ey) {xie y -\- zx2ey) 

Substituting these expressions in eq. ( |B.lD and eq. ( [B.2D one obtains 

Ml{h,l2,k) 



C^-{tl+Zt2). 



^^^'■g^il + z') [A,A,{tl + zHl) + BiBq{tl-zhl)] , 



zk^ 



Nr 



(B.6) 
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and 

MUh,l2, k)^^g'^{l-z)[z^ + il- zf] [AiA,itl + zhl) + B,B^{t\ - zHl)\ . 

t ^ 

(B.7) 

where the colhnear singularity has been isolated. Note that the result is proportional 
to the Born matrix element in eq. (|3.1| ) with X2 replaced by ZX2, indicating that the 
momentum fraction of the parton entering the hard scattering has been reduced by a 
factor z after the emission of a collinear gluon. 
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